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ABSTRACT 


We  have  investigated  the  characteristics  of  P-coda  and  Lg 
measurements  at  the  NORSAR)  and  Graefenburg,  arrays*  for  presumed 

underground  nuclear  explosions  in  the  Semipalatinsk  region  of  the 

T  u 

Soviet  Union.  Our  main  objectives  -in  this  study>were  to  investi¬ 
gate  the  effects  of  the  propagation  paths  in  western  Russia  on 
the  narrowband  and  broadband  recordings  of  Lg  at  teleseismic 
distances  and  to  study  the  relative  P-coda  and  Lg  amplitudes 
recorded  at  these  two  arrays  for  the  largest  ffmh,  ■>  6.0)  Shagan 
River  explosions.  '  ,  ,  .  -  ;  - 

y  ( •  •-  •  ) 

\ 

—  ^Comparison  of  broadband  recordings  of  teleseismic  Lg  at 
Graef  enburq^fA  =  42Q  V1  with  narrowband  NORSAR  /,  ( A=38°)  ‘  and  filtered 
Graefenburg  recordings  of  Lg  from  Shagan  River  events  reveals 
that  Lg  is  more  obvious,  relative  to  the  preceding  P-coda,  on 
broadband  seismograms  than  on  high-frequency  seismograms.  Broad¬ 
band  recordings  of  Lg  at  Graefenburg  are  about  0.5  log  units 
stronger  in  the  0.2  -  1.0  Hz  band  than  in  the  0.6  -  3.0  Hz  range 
although  noise  is  also  correspondingly  higher. 


The  early  P-coda  at  NORSAR  is  stronger,  relative  to  Lg,  than 
that  at  Graefenburg.  Also,  the  coda-envelope  shapes  are  quite 
different  for  the  two  arrays^^the  NORSAR  coda  flattens  between 
about  200  to  340  seconds  after  P,  whereas  the  Graefenburg  coda 
decays  monitonically  from  the  P-wave  out  to  the  Sn  and  Lg  arrival 
times.  The  standard  deviations  of  Lg  and  S-coda  phases  at  NORSAR 
are  higher  by  about  0.05  to  0.06  log-rms  amplitude  units  than 
those  of  P-coda  phases,  which  is  probably  related  to  differential 
site  and  propagation  effects  on  P,  S  and  Lg  waves.  The  most 
stable  part  of  the  P-coda  before  Lg  is  the  flat  part  of  the  coda, 


where  the  spatial  standard  deviation  in  log-rms  amplitude  in  5 
second  windows  drops  to  a  minimum  of  less  than  0.1  units.  The 
flat  coda  envelopes  from  about  200  to  340  seconds  after  P, 
recorded  at  NORSAR  for  Semipalatinsk  explosions,  are  consistent, 
on  the  basis  of  timing,  to  Lg  -  P  forward  scattering  from  lateral 
heterogeneities  in  the  Ural  Mountains. 

Measurements  on  Graefenburg  data  of  vertical  and  transverse 
log-rms  amplitudes  in  the  P-coda  and  Lg  phase  in  two  frequency 
bands  correlate  well  with  network  mb.  These  measurements  indi¬ 
cate  that  Lg  and  P-coda  amplitudes  do  not  appear  to  be  strongly 
affected  by  the  tectonic  component  as  revealed  by  surface-wave 
studies.  However,  the  slopes  of  the  regression  of  Lg  amplitude 
versus  mb  were  small  ranging  from  0.6  to  0.8.  The  slopes  of  the 
P-coda  regressions  were  higher,  ranging  from  0.9  to  over  1.0. 
This  supports  the  idea  that  early  P-coda  waves  are  generated  by  P 
scattering.  The  principal  cause  of  the  small  slopes  is  that  the 
variance  of  the  Lg  and  P-coda  amplitudes  for  the  largest  Shagan 
River  events  are  significantly  smaller  than  that  of  the  corres¬ 
ponding  network-averaged  mbs.  We  argue  that  P-coda  and  Lg  ampli¬ 
tudes  are  more  robust  estimators  of  relative  yield  than  mb,  and 
that  the  yield  range  for  these  events,  based  on  Lg  and  P-coda 
measurements  at  NORSAR  and  Graefenburg,  is  smaller  by  approxi¬ 
mately  a  factor  of  2  than  the  yield  range  indicated  by  mb.  This 
is  consistent  with  the  conjecture  of  Sykes  and  Cifuentes  (1983) , 
based  on  surface-wave  studies,  that  the  largest  Shagan  explo¬ 
sions,  with  mb  greater  than  6.0,  have  nearly  the  same  yield. 
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Since  mid- 1978,  the  body-wave  and  surface-wave  magnitudes  of 
the  largest  Soviet  nuclear-weapons  tests  at  the  Shagan  River  test 
site  have  increased  by  nearly  a  factor  of  two  (Alewine  and  Bache, 
1983).  Typically,  the  body-wave  magnitudes,  m^,  for  Shagan  River 
explosions  fired  since  1978  Have  been  in  the  range  of  5.8  to  6.2. 
Because  explosions  with  thes..  magnitudes  detonated  at  the  Nevada 
Test  Site  (NTS)  have  yields  well  in  excess  of  150  kt,  the  large 
m^s  of  the  recent  Shagan  River  explosions  have  raised  concerns 
about  possible  violations  of  the  Threshold  Test  Ban  Treaty 
(TTBT) . 

The  accuracy  of  absolute  yield  estimates  of  underground 
nuclear  explosions,  using  seismic  methods,  is  reduced  by  a  number 
of  factors,  the  most  important  one  being  the  uncertainty  in  the 
attenuation  bias  between  NTS  and  the  Shagan  River  test  site. 
Several  direct  and  indirect  methods  have  been  applied  to  deter¬ 
mine  the  m^  bias  so  that  the  m^-yield  calibration,  determined  for 
NTS,  can  be  directly  used  to  estimate  yields  of  explosions  at 
Shagan  River.  However,  the  large  range  of  bias  estimates  (0.2  to 
0.4  m^)  makes  it  difficult  to  determine  with  the  yields  ot  the 
largest  Shagan  explosion  with  a  large  degree  of  certainty. 

An  alternative  approach  is  to  utilize  measurements  which  are 
not  as  much  affected  by  attenuation  bias  as  m^.  One  such 
measurement  is  the  surface-wave  magnitude,  or  Mg .  Sykes  and 
Cifuentes  (1984)  estimated  the  yields  of  the  largest  Shagan  River 
explosions  with  low  tectonic  release  with  Mg  measurements  and 
concluded  that  the  explosions  with  m^s  ranging  from  6.0  to  6.2 
had  essentially  the  same  yield  but  less  than  150  kt.  They  argued 


that  the  large  variance  in  results  from  geologic  variations  in 
the  source  region  and  measurement  uncertainties. 

The  problem  with  Mg  yield  estimation  is  that  surface  waves 
can  be  affected  by  non-isotropic,  tectonic  component.  Sykes  and 
Cifuentes  (1984)  avoided  this  problem  by  selecting  events  for 
calibration  and  yield  estimation  which  had  low  tectonic  compo¬ 
nent.  Given  et  al  (1983)  and  Given  and  Mellman  (1984)  have 
applied  source  inversion  techniques  to  th*5  Shagan  River  explo¬ 
sions  in  order  to  obtain  a  magnitude  correction  for  the  tectonic 
component.  However,  this  method  is  still  being  tested. 

The  Lg  phase  has  proven  to  be  useful  for  yield  estimation 
because  it  combines  the  advantages  of  surface  waves  along  with 
the  fact  that  Lg  magnitudes  appear  to  be  insensitive  to  tectonic 
component  (Alexander,  1984).  Nuttli  (1984)  has  developed  an 
absolute  Lg  magnitude  method  in  which  Lg  magnitude  attenuation 
corrections  are  determined  from  Lg-coda  Q  estimates.  Although 
Nuttli  has  had  great  success  with  applying  this  method  to  events 
of  known  yield,  other  investigators  have  not  been  as  successful 
(Rondout  Associates,  1984).  The  problem  seems  to  be  large  uncer¬ 
tainties  in  the  estimation  of  Lg-coda  Q  from  the  time  decay  of 
frequencies  in  the  coda.  Nuttli  (1984)  has  suggested  that  the 
possible  causes  of  these  uncertainties  are  interference  of  funda¬ 
mental  mode  Rayleigh  waves  and  noise  with  the  coda  which  can 
cause  inaccurate  frequency  estimates  at  large  epoch  times  into 
the  coda. 

Nuttli  (1984)  has  applied  his  Lg  magnitude  yield  estimation 
technique  to  regional  recordings  of  Shagan  River  explosions, 
including  some  of  the  largest  explosions.  His  yield  estimates 
for  the  largest  explosions  (mfc>  between  6.0  and  6.2)  range  from 
120  to  229  kt.  Also,  by  comparing  Lg  magnitudes  with  network- 
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averaged  mb  for  events  at  NTS  and  Shagan  River  and  assuming 
absolute  Lg  magnitudes  to  be  insensitive  attenuated  bias,  Nuttli 
(1984)  obtains  a  preliminary  mb  bias  between  NTS  and  Shagan  River 
of  about  0.41  magnitude  units. 

The  network  bias  problem  could  in  principle  be  eliminated  if 
the  TTBT  were  ratified,  because  the  protocol  of  the  treaty  calls 
for  two  calibration  shots  of  known  yield  to  be  fired  within  each 
distinct  geologic  region  at  each  test  site  in  the  U.S.  and  the 
USSR.  Assuming  that  the  yields  of  the  calibration  explosions  are 
accurate,  the  mb  attenuation  bias  could  be  directly  determined. 
However,  there  would  still  remain  the  question  of  the  precision 
of  mb  estimates  for  relative  yield  estimation.  Sykes  and 
Cifuentes  (1984)  have  suggested  an  mb  variance  as  great  as  0.2 
magnitude  units  for  largest  Shagan  River  explosions  of  the  same 
yield  which  could  result  in  close  to  a  factor  of  two  uncertainty 
of  yield  estimation. 

There  have  been  a  number  of  studies  which  have  indicated 
that  measurements  made  in  the  P-coda  are  at  least  as  stable  as 
network-averaged  mb  estimates  and  significantly  better  than 
single-station  mb  (Baumgardt,  1983;  Bullitt  and  Cormier,  1984; 
Baumgardt,  1984;  Gupta  et  al,  1984).  Similar  results  have  been 
obtained  for  Lg  and  Lg-coda  wave  measurements  (Ringdal,  1983; 
Alexander,  1984).  These  results  suggest  that  the  yield  estima¬ 
tion  precision  may  be  improved  by  using  magnitudes  derived  from 
long  averaging  windows  in  the  P-coda  and  on  the  Lg  phase. 
Apparently,  such  measurements  average  out  the  effects  of  focusing 
and  defocusing  which  cause  the  large  scatter  in  P-wave  magni¬ 
tudes  . 

In  this  report,  we  present  the  results  of  a  study  of  the 
precision  of  using  P-coda  and  Lg  measurements  for  relative  yield 


estimation.  In  particular,  we  focus  on  the  P-codas  and  Lg-waves 
recorded  at  two  seismic  arrays,  NORSAR  and  Graefenburg,  from  the 
largest  underground  explosions  detonated  at  Shagan  River.  Table 
1  gives  a  list  of  the  events,  their  m^s  as  determined  by  Sykes 
and  Cifuentes  (1984)  and  NEIS,  and  the  NORSAR  and  Graefenburg 
data  which  were  available  for  this  study.  In  the  first  part  of 
the  report  we  discuss  the  broad  characteristics  of  the  long  P- 
codas  recorded  at  these  two  arrays  for  the  largest  Shagan  River 
explosions.  In  this  section,  we  suggest  that  Lg  forward  scat¬ 
tering  from  heterogeneities  along  the  path  between  the  source  and 
receiver  may  be  important  for  generating  P-coda  waves.  In  the 
second  part  of  this  report,  we  investigate  the  conjecture  of 
Sykes  and  Cifuentes  (1984)  that  the  Shagan  River  explosions  with 
mbs  greater  than  6.0  are  of  approximately  the  same  yield,  by 
analyzing  the  relative  P-coda  and  Lg  amplitudes  for  these  events. 
Our  overall  conclusions  from  this  study  are  that  single-site  P- 
coda  and  Lg  measurements  are  as  precise  as  network  Ms  measure¬ 
ments  for  relative  yield  estimation  and  indicate  that  the  largest 
Shagan  River  explosions  are  closer  in  yield  than  indicated  by 
their  network  m^  estimates. 


TABLE  1 

SHAGAN  RIVER  DATA  USED  IN  THIS  STUDY 


Date 

* 

mb 

NORSAR 

Data 

GRAEFENBURG 

Data 

15  Sep  1978 

5.963 

X 

X 

4  Nov  1973 

5.576 

X 

X 

29  Nov  1978 

5.996 

X 

23  Jun  1979 

6.215 

X 

X 

7  Jul  1979 

5.839 

X 

X 

4  Aug  1979 

6.161 

X 

X 

18  Aug  1979 

6.170 

X 

28  Oct  1979 

5.990 

X 

2  Dec  1979 

5.998 

X 

23  Dec  1979 

6.170 

X 

12  Jun  1980 

5.6 

X 

29  Jun  1980 

5.707 

X 

12  Oct  1980 

5.918 

X 

X 

14  Dec  1980 

5.953 

X 

13  Sep  1981 

6.064 

X 

29  Mar  1981 

5.6 

X 

22  Apr  1981 

5.954 

X 

25  Apr  1982 

6.1 

X 

5  Dec  1982 

6.1 

X 

26  Dec  1982 

5.7 

X 

12  Jun  1983 

6.1 

X 

6  Oct  1983 

6.0 

X 

1 

26  Oct  1983 

6.1 

X 

19  Feb  1984 

5.9 

X 

15  Apr  1984 

5.7 

X 

25  Apr  1984 

6.0 

X 

26  May  1984 

6.1 

X 

*  Magnitude  estimates  to  two  3-decimal  places  are  from  Sykes  and 
Cifuentes  (1983).  All  others  are  from  NEIS. 


2.0  TELESEISMIC  P-CODA 
AND  Lg-WAVES  IN  EURASIA 


2.1  INTRODUCTION 

In  seismological  studies  of  earthquakes  and  quarry  blasts  at 
local  and  regional  distances,  the  short-period  P  and  Pn  coda,  Lg, 
and  Lg  coda  have  played  important  roles.  Theoretical  and  empir¬ 
ical  attenuation  ..adies  by  means  of  spectral  analysis  of  P  and 
Lg  codas  have  provided  insights  into  the  relative  importance  of 
intrinsic  and  scattering  attenuation  of  Lg  and  lateral  variations 
of  attenuation  in  the  earth's  crust  and  lithosphere.  (Herrmann, 
1980;  Singh  and  Herrmann,  1983;  Herrmann  and  Wang,  1983.)  More¬ 
over,  P-coda  duration  magnitude  and  m^  (Lg)  magnitudes  have  been 
routinely  determined  for  local  and  regional  seismic  events  and 
have  long  been  known  to  provide  very  stable  magnitude  estimates. 
(Herrmann,  1975;  Bakun  and  Lindh,  1977;  Suteau  and  Whitcomb, 
1979;  Shapira,  1981.) 

In  the  area  of  seismological  monitoring  of  nuclear  test  ban 
treaties,  yield  estimation  for  verification  of  the  threshold 
treaty  is  based  almost  entirely  on  mb  estimates  at  teleseismic 
distances.  However,  recent  studies,  discussed  above,  have 
revealed  that  mb  measures  are  subject  to  a  number  of  biases,  one 
of  the  most  important  being  focusing  and  defocusing  of  short- 
period  P-waves  by  lateral  heterogeneities  in  the  earth.  Tele¬ 
seismic  P-coda  and  regional  Lg  magnitudes  have  proven  to  be 
surprisingly  stable  in  spite  of  the  complexity  of  high-frequency 
scattered  and  higher-mode  waves. 


In  this  section,  we  discuss  the  characteristics  of  long-term 
codas  (P,  P-coda,  Lg,  Lg-coda)  recorded  at  the  two  arrays, 
Graefenburg  and  NORSAR,  from  the  largest  Shagan  River  explosions. 
NORSAR  and  Graefenburg  are  at  distances  of  38°  and  42°,  respec¬ 
tively,  and  thus,  are  well  into  the  teleseismic  distance  range. 
In  this  section  we  show  that  Lg  is  well  recorded  at  these  two 
arrays.  Further,  we  suggest  that  forward  scattering  of  Lg-to-P 
from  heterogeneities  between  the  source  and  receiver  may  generate 
some  of  the  pre-Lg  P-coda  waves.  The  results  are  important  in 
regard  to  our  understanding  of  the  origin  of  the  P-coda,  the 
nature  of  scattering  attenuation  of  Lg ,  and  the  utility  of  P-coda 
and  Lg  measurements  for  yield  estimation. 

2.2  DATA  AND  ANALYSIS  PROCEDURES 

Figure  1  shows  the  location  of  the  Semipalatinsk  test  area 
and  the  two  arrays,  Graefenburg  and  NORSAR.  The  propagation 
paths  from  the  Semipalatinsk  test  site  to  both  of  these  arrays 
are  very  similar,  although  there  may  be  some  near  receiver 
differences . 

Propagation  paths  in  western  Russia  are  primarily  shield 
type  and  Lg-wave  propagation  across  this  region  has  been  found  to 
be  very  efficient  (North,  1978;  Gupta  et  al,  1980;  Nuttli,  1981). 
North  (1978)  found  that  the  observed  attenuation  of  Lg-waves  at 
Scandinavian  stations  from  sources  distributed  throughout  western 
Russia  was  roughly  consistent  with  that  determined  by  Nuttli 
(1973)  for  eastern  United  States,  orr=0.07  deg-1  *-tli  (1981) 
observed  a  somewhat  higher  value  of  y^O.133  de  (  or  similar 
paths  which,  nevertheless,  is  still  quite  low.  Noi  t-h  (1978)  also 
concluded  that  the  Ural  Mountains  and  Baltic  Sea  do  not  block  or 
strongly  attenuate  Lg-waves  which  cross  these  structures  from 
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FIGURE  1  MAP  OF  WESTERN  RUSSIA 


events  in  western  Russia  to  stations  in  Scandinavia.  Thus,  the 
low  attenuation  ot  Lg  along  the  stable  Baltic  shield  paths  ac¬ 
counts  for  the  teleseismic  observation  of  Lg  at  the  Scandinavian 
stations . 

The  configuration  of  the  NORSAR  array  is  illustrated  in 
Figure  2.  The  array  originally  consisted  of  22  subarrays,  with 
each  subarray  consisting  of  six  short-period,  vertical  instru¬ 
ments.  The  horizontal  diameter  was  100km.  Since  October  1, 
1976,  the  array  has  consisted  of  seven  subarrays  spread  over  a  50 
km  aperture  and  shown  in  Figure  2  by  the  filled  circles. 

The  Graef enburg-array  instrument  locations  are  shown  in 
Figure  3.  In  the  present  study  we  have  concentrated  on  data 
recorded  at  the  Al  three-component  sensor,  which  is  also  the 
location  of  the  GRFO  SRO  station.  The  unique  feature  of  the 
Graefenburg  array  data  is  that  it  is  recorded  broadband,  with  an 
instrument  response  flat  to  velocity  from  about  20  seconds  period 
to  5  Hz  (Harjes  and  Seidl,  1978).  In  contrast,  the  NORSAR 
sensors  record  data  in  the  traditional  short-period  and  long- 
period  bands.  The  mid-period  band,  from  about  1.6  to  5  seconds 
is  not  sampled  at  NORSAR,  but  is  passed  by  the  broadband 
instruments  at  Graefenburg. 

Table  1  lists  the  events  which  were  examined  in  this  study, 
their  magnitudes  as  determined  by  Sykes  and  Cifuentes  (1983)  and 
NEIS,  and  the  available  data  from  NORSAR  and  Graefenburg.  Table 
1  shows  that  we  have  data  for  several  events  that  were  recorded 
by  both  NORSAR  and  Graefenburg.  Comparison  of  seismograms 
recorded  at  the  two  arrays  for  common  events  allows  us  to  examine 
the  effect  of  the  bandpasses  on  the  coda  and  Lg  recordings  and  to 
look  for  subtle  differences  between  the  Shagan-NORSAR  and  Shagan- 
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FIGURE  2  DEPLOYMENT  OF  THE  ORIGINAL  22  NORSAR  SUBARRAYS  (FROM  1971  TO 

OCTOBER  1,  1976).  THE  SUBSET  OR  7  SUBARRAYS  IN  OPERATION  SINCE 
OCTORF.R  1,  1976  ,  APE  MARKED  AS  FILLED  CIRCLES. 
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LOCATION  OF  GRAEFENBURG  ARRAY  STATIONS  (FROM 
HARJES  AND  SEIDL ,  19  78  ) 
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FIGURE  4  TIME  WINDOWS  USED  FOR  NOISE  AND  P-CODA  MAGNITUDE  MEASUREMENTS 
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In  our  analysis  of  long-term  P-codas,  we  have  found  it  very 
useful  to  plot,  as  coda  envelopes,  the  log-rms  amplitudes  as  a 
function  of  time  in  adjacent  5  second  windows.  For  a  single 
channel,  the  signal  level  in  each  window  is: 


Sjk  =  lQ9  Ajk 


(3) 


for  the  k'th  channel.  The  average  across  an  array  is: 
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and  the  standard  deviation  is 


SD?  ■ 


(5] 


where  L  is  the  number  of  channels  in  the  array.  For  a  single 
channel,  we  compare  Sj^  as  a  function  of  time,  corresponding  to 
the  start  of  each  window,  with  the  average  noise  level,  R^.  in 
the  multichannel  case,  we  compare  the  Sj  values  against  the  noise 
levels  averaged  across  the  array: 


Notice  that  in  (6)  we  have  taken  the  channel  average  of  log-rms 
amplitude  before  taking  the  average  over  all  the  windows,  instead 
of  averaging  the  window  average  in  (2)  across  all  channels  in  the 
array.  The  corresponding  noise  standard  deviation  is 


-  E 


(log  A. 


Thus,  in  comparing  Sj  with  N,  ,SDjS  with  SDn  and  Sjk  with  Rk,  we 
are  comparing  each  5  second  windowed,  log-rms  P-coda  amplitude  as 
a  function  of  time  with  the  average  of  the  log-rms  noise  ampli¬ 
tude,  in  5  second  windows,  over  the  50  seconds  of  noise  ahead  of 


We  determine  average  P-coda  and  Lg  magnitude  in  the  same  way 
that  we  compute  the  average  noise  levels  for  a  single  channel  in 
(2) .  Thus . 


where  Aj^  is  an  rms  amplitude  measurement  in  the  j ' th  window  and 
the  k'th  channel,  and  Mg  is  the  number  of  windows. 

In  past  stjdies,  (Baumgardt,  1983,  1984),  we  have  also 
examined  array-average  coda  and  Lg  log-rms  amplitudes  determined 
in  the  same  manner  as  the  noise  in  (6).  However,  Ringdal  (1983) 
and  Baumgardt  (1984)  showed  that  array-average  measurements  were 
almost  perfectly  correlated  with  single-channel  measurements. 
Thus,  in  this  study,  we  have  only  made  single-channel  measure¬ 
ments  . 

2.3  COMPARISON  OF  P-CODA  AND  Lg  RECORDINGS  AT  NORSAR  AND 
GRAEFENBURG 

Figure  5  shows  traces  of  three  large  Semipalatinsk  explo¬ 
sions  recorded  at  one  of  the  NORSAR  sensors.  The  plots  are 
highly  compressed,  with  each  plot  covering  a  time  span  of  20 
minutes,  and  the  tick  marks  on  the  time  scale  denote  10  seconds. 
The  plots  reveal  the  emergent  nature  of  the  teleseismic  Lg  phase 
recorded  at  NORSAR.  It  would  be  very  difficult  to  find  Lg  on 
these  traces  if  they  weren't  plotted  in  this  highly  compressed 
form.  Also,  note  the  variability  of  the  relative  amplitude  of 
the  P  and  PP  phases  for  these  events,  which  are  only  30  to  40  km 
apart.  This  is  an  example  of  the  effect  of  ray-parameter  depen¬ 
dent  focusing  and  defocusing  in  the  earth  which  causes  signifi¬ 
cant  variations  of  P-wave  magnitudes. 

A  Graefenburg  recording  of  a  Shagan  River  explosion  is  shown 
in  Figure  6.  The  top  trace  is  the  original  broadband  trace  and 
the  bottom  trace  is  the  broadband  seismogram  passed  through  a 
high-pass  filter.  The  filter  is  a  6-pole  recursive  filter  with  a 
low-frequency  cutoff  of  0.6  Hz,  which  is  close  to  that  of  the 
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FIGURE  b  BROADBAND  (TOR)  AND  FILTERED  ^BOTTOM)  TRACES  OF  THE  4  AUGUST  1979 
SHAGAN  RIVER  EVENT  RECORDED  AT  THE  GRAEFENBURG  ARRAY 


NORSAR  instrument  filter.  Notice  that  the  Lg  wave  on  the  broad¬ 
band  trace  is  very  strong  and  appears  to  be  a  distinct  phase  as 
opposed  to  a  slight  increase  in  the  coda  level,  as  in  Figure  6. 
However,  Lg  is  less  clear  on  the  bottom  trace  which  is  filtered 
to  resemble  the  NORSAR  instrument  bandpass. 

Figure  7  compares  seismograms  written  at  NORSAR  and  Graefen- 
burg  by  an  explosion  at  Shagan  River.  Notice  that  the  Lg  wave  is 
more  obvious  on  the  broadband  Graefenburg  trace  than  on  the  high- 
frequency  NORSAR  trace.  The  third  trace  at  the  bottom  of  Figure 
7  was  made  by  passing  the  broadband  Graefenburg  trace  through  the 
high-pass  filter  which  simulates  the  high-frequency  NORSAR  re¬ 
sponse.  In  the  high-frequency  band  at  Graefenburg,  the  Lg  phase 
is  weaker  relative  to  the  preceding  coda  than  in  the  broadband 
recording.  Also,  comparing  the  top  and  bottom  traces  reveals 
that  the  coda  ahead  of  the  Lg  phases  appears  to  be  more  intense 
at  NORSAR  than  at  Graefenburg.  Also,  note  that  although  the  Lg 
phase  is  stronger  in  the  intermediate  band  than  in  the  high- 
frequency  band,  so  also  is  the  noise  level  ahead  of  the  first 
arrival  P. 

Another  representation  of  explosion  codas  recorded  at  the 
NORSAR  and  Graefenburg  arrays  is  given  in  Figure  8,  where  we  have 
plotted  single-channel  log-rms  amplitudes,  computed  in  5  second 
windows  as  a  function  of  time,  as  discussed  above  (Equation  3). 
Each  seismogram  was  passed  through  a  0.6  -  3.0  Hz  filter  prior  to 
computing  the  rms  amplitudes.  The  horizontal  lines  on  each  plot 
denote  the  average  rms  amplitude  in  the  50  seconds  of  noise 
background  ahead  of  the  first-arrival  P.  Thus,  we  are  comparing 
5  second  coda  levels,  determined  by  equation  3,  with  the  average 
5  second  noise  level,  determined  by  equation  2,  in  the  50  seconds 
of  noise  ahead  of  P.  These  plots  show  that  the  Lg  amplitudes, 
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FIGURE  8  NORSAR  (TOP)  AND  GRAEFENBURG  (BOTTOM)  CODA  ENVELOPES  OF 
THE  4  AUGUST  1979  AND  22  APRIL  1981  SHAGAN  RIVER  EVENTS 


relative  to  the  noise,  are  about  the  same  at  NORSAR  and  Graefen- 
burg,  which  is  consistent  with  the  observation  in  Figure  7.  The 
signal-to-noi se  ratios  for  Lg  at  each  array  are  about  12  to  14  dB 
in  the  0.6  -  3.0  Hz  band.  Also,  note  that  the  pre-Lg  coda  level 
appears  slightly  greater,  relative  to  noise,  at  NORSAR  than  at 
Graef enburg . 

Figure  9  shows  coda  plots  for  a  Shagan  River  explosion 
recorded  at  Graefenburg  in  two  frequency  bands:  0.2  -  1.0  Hz  at 
top  and  0.6  -  3.0  Hz  at  bottom.  Comparison  of  these  two  plots 
shows  that  the  Lg  level  is  clearly  higher  relative  to  the  coda 
level  in  the  mid-period  band  than  in  the  short-period  band. 
However,  note  that  the  noise  level  is  also  higher  in  the  mid¬ 
period  band  than  in  the  short-period  band,  which  results  in  the 
Lg  signal-to-noi se  ratio  being  about  the  same  in  the  two  bands. 

These  results  illustrate  the  advantage  and  disadvantages  of 
broadband  recording  of  regional  seismic  phases.  Figures  6,  7, 
and  8  show  that  Lg  excitation  is  stronger  in  the  mid-period  range 
than  in  the  high-frequency  band,  a  fact  also  pointed  out  by 
(Noponen  et  al,  1979).  However,  the  noise  level  is  also  higher. 
Thus,  broadband  signal-to-noi se  ratios  are  not  significantly 
higher  than  those  in  the  high-frequency  band  although  broadband 
Lgs  do  stand  out  more  clearly  than  high-frequency  Lgs.  Thus,  Lg 
magnitudes  may  be  measured  more  easily  and  reliably  on  broadband 
traces  than  on  high-frequency  traces.  However,  another  possible 
problem  with  low-frequency  Lg  magnitudes  is  that  they  may  be  more 
affected  by  tectonic  component  than  high-frequency  Lg  magnitudes. 
We  shall  consider  this  point  again  in  Section  3.0. 
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FIGURE  9  MID-BAND  (TOP)  AND  HIGH-BAND  (BOTTOM)  FILTERED  CODA  ENVELOPES 
FROM  23  JUNE  1979  SHAGAN  RIVER  EVENT  RECORDED  AT  GRAEFENBURG 


2.4 


SPATIAL  VARIATION  ACROSS  NORSAR  OF  P-CODA  AND  Lg-WAVE 
AMPLITUDES 


Figure  10  shows  coda  plots  for  16  Semipalatinsk  explosions 
recorded  at  the  01A01  channel  at  NORSAR.  Clearly,  the  plots 
reveal  significant  scatter  caused  by  local  scattering  beneath  the 
receiver.  These  effects  are  essentially  random,  and  the  variance 
can  be  reduced  by  averaging  the  rms  amplitudes  across  all  avail¬ 
able  channels  (equation  4)  .  The  result  of  this  processing  is 
shown  in  Figure  11,  where  the  time  variance  has  been  reduced 
enough  such  that  a  structure  can  be  discerned.  All  codas  show  a 
noticeable  flattening  starting  at  about  200  seconds  and  extending 
to  about  340  seconds.  A  second  flattening  starts  at  about  480 
seconds  and  extends  to  about  the  Lg  time  of  about  745  seconds. 

Figure  12  shows  a  plot  of  the  standard  deviations  in  the 
log-rms  amplitudes  across  the  array  as  a  function  of  time  (equa¬ 
tion  5) .  All  the  codas  exhibit  a  drop  in  standard  deviation  from 
about  0.2  to  0.25  for  the  P  and  PP  waves  to  0.05  to  0.1  in  the 
early  coda.  Notice  that  the  noise  standard  deviation  in  the  log- 
rms  levels  is  about  0.12.  The  standard  deviations  begin  rising 
at  about  320  seconds,  peak  at  about  0.15  at  800  seconds,  and  tail 
off  thereafter. 

Figure  13  shows  the  log-rms  amplitudes,  averaged  across  the 
array,  and  their  standard  deviations  plotted  with  time  for  three 
Shagan  River  events.  The  large  standard  deviations  associated 
with  the  P  and  PP  arrivals  are  clearly  evident  and  are  probably  a 
result  of  focusing  and  defocusing  effects  from  laterally  hetero¬ 
geneous  structure  beneath  the  array.  Another  phase  of  interest, 
which  arrives  about  60  seconds  after  P,  is  labelled  P420 
Figure  13.  This  phase  is  very  coherent  across  the  array  and 
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appears  on  almost  ail  or  tne  array-averaged  coda  envelopes  for 
Semipalatinsk  explosions  (see  Figure  11)  .  Like  P  and  PP,  its 
standard  deviation  around  the  array  is  higher  than  the  other  coda 
waves. 

We  suggest  that  this  phase  is  a  wide-angle,  post-critical 
reflected  or  diffracted  phase  from  the  top  of  the  420  kilometer 
discontinuity,  which  is  the  reason  for  our  calling  the  phase 
p420*  Rin9  and  Calcagnile's  (1976)  mantle  model  KCA  for  western 
Russia,  whose  travel  time  curve  is  shown  in  Figure  14,  predicts 
such  wide-angle  reflections  at  distances  near  38°.  King  and 
Calcagnile  (1976) ,  in  fact,  observed  wide-angle  reflections  out 
to  about  32°,  but  extrapolated  the  be  branch  of  their  travel  time 
curve  out  to  near  38°,  as  shown  in  Figure  14.  Based  on  this 
extrapolation,  the  reflected  phase  at  38°  should  arrive  about  50 
seconds  after  P  and  have  a  slowness  of  about  12  seconds/degree. 
Using  the  frequency-wave  number  method,  we  measure  a  slowness  of 
12.6  seconds/degree  for  the  phase  P42o*  However,  as  we  mentioned 
above,  the  phase  arrives  about  10  seconds  later  than  predicted  by 
the  model  KCA,  or  60  seconds  after  P.  Another  possibility  is 
that  the  phase  may  be  associated  with  a  group  of  deterministi¬ 
cally  scattered  phases  from  the  420  discontinuity  referred  to  by 
Menke  and  Richards  (1982)  as  "whispering  gallery  phases,"  which 
may  be  expected  to  arrive  later  than  the  topside,  wide-angle 
reflected  phase. 

Aside  from  these  phases,  most  of  the  rest  of  the  coda  out  to 
the  Lg  phase  seems  to  be  non-determini stically  scattered  phases. 
As  has  been  observed  in  earlier  studies  (Nersesov  et  al,  1975; 
Baumgardt,  1983;  Bullitt  and  Cormier,  1984;  Baumgardt,  1984), 
high-frequency  scattered  waves  in  the  P-coda  have  lower  spatial 
variance  than  deterministic  phases.  In  the  early  coda  between  P 


and  PP  and  in  the  flat  part  of  the  coda  between  about  160  -  200 
seconds  to  340  seconds,  in  Figure  13,  the  standard  deviations 
reach  their  lowest  levels  of  about  0.08  to  0.1.  Starting  at 
about  the  expected  arrival  time  of  S-waves,  the  standard 
deviations  begin  to  increase.  It  is  interesting  that  there  is  no 
evidence  of  increased  amplitudes  as  a  result  of  S-wave  arrivals; 
in  fact,  the  coda  level  actually  starts  to  decrease  at  that 
time.  Thus,  increased  standard  deviation  may  be  the  strongest 
indication  of  the  arrival  of  shear  -  wave  coda  waves.  The 
standard  deviation  peaks  at  about  0.15  at  the  Lg  arrival  time  and 
then  falls  off  again  in  the  Lg-coda. 


What  is  the  cause  of  the  high  spatial  standard  deviation  in 
the  Lg-wave?  Given  the  observed  stability  of  Lg-wave  magnitudes, 
it  was  expected  that  the  standard  deviations  would  be  quite  low. 
In  fact,  Ringdal  (1983)  obtained  lower  standard  deviations  (0.06 
-  0.08  log-rms  units)  than  those  in  Figures  12  and  13.  The 
apparent  reason  for  this  is  that  Ringdal' s  measurements  on  Lg 
were  made  in  2  minute  windows,  whereas  those  in  Figures  12  and  13 
were  made  in  5  second  windows.  However,  the  fact  still  remains 
that  the  short-window  Lg  standard  deviations  are  at  least  0.05 
units  higher  than  those  of  the  earlier  P-coda  and  later  Lg-coda. 


Baumgardt  (1983)  suggested  that  the  higher  standard  devia¬ 
tions  of  Lg  may  be  explained  by  the  fact  that  the  array  averages 
were  determined  after  first  lining  up  the  traces  on  the  first 
arrival  P-waves.  Thus,  because  Lg  propagates  more  slowly  than  P 
(Lg  apparent  velocity  would  be  3.5  to  4.5  km/sec  compared  with 
the  expected  P-wave  apparent  velocity  of  about  13.21  km/sec),  the 
Lg-waves  would  be  misaligned  on  what  are  essentially  incoherent 
P-wave  beams  in  Figures  12  and  13.  To  check  this,  we  compare,  in 
Figure  15,  log-rms  standard  deviation  plots  for  a  Semipalatinsk 


event  with  the  codas  aligned  to  a  P-wave  (Figure  15a)  and  an  Lg- 
wave  (Figure  15b)  velocity.  The  Lg  velocity  was  assumed  to  be 
3.5  km/sec.  The  comparison  reveals  a  higher  noise  standard 
deviation  for  the  Lg  velocity  average  because  the  P  wave,  which 
is  not  lined  up  correctly,  contaminates  the  background  noise 
estimate.  However  the  coda  shapes  for  the  P  and  Lg  beams  are 
almost  identical.  Even  when  the  codas  are  lined  up  consistent 
with  an  Lg  velocity,  the  Lg  log-rms  standard  deviation  is  still 
0.03  to  0.05  units  higher  than  that  of  the  earlier  P  coda. 

It  appears  then  that  higher  standard  deviation  is  an  intrin¬ 
sic  characteristic  of  the  Lg  phase  and  its  coda.  Given  the  low 
coherence  of  Lg-waves  over  intersensor  separation  as  small  as  2 
to  5  km  (Mrazek  et  al,  1980;  Mykkeltveit  et  al,  1980;  Der  efc  al, 
1984)  it  is  probably  not  unsurprising  that  the  Lg  standard  devia¬ 
tions  across  the  entire  NORSAR  array,  with  aperture  of  about  50 
km,  do  not  depend  strongly  on  the  beam  velocity.  Also  the  higher 
standard  deviations  of  Lg  in  Figures  12  and  13  last  for  several 
seconds  after  the  3.5  km/sec  group  velocity  time,  which  is  much 
longer  than  the  time  into  the  Lg  coda  where  any  spatial  coherence 
has  been  observed  (Der  et  al,  1984). 

We  argue  that  the  high  standard  deviation  in  the  Lg  part  of 
the  coda  is  caused  by  lateral  variations  local  structure  beneath 
the  sensors  in  the  NORSAR  array.  The  fact  that  Lg  excitation  is 
highly  sensitive  to  geologic  variations  and  crustal  structure  has 
been  well  documented  (Barker  et  al,  1981).  It  is  also  interest¬ 
ing  that  the  standard  deviations  in  the  coda  begin  to  increase  at 
about  the  expected  arrival  times  of  S  waves,  even  though  no 
strong  S-waves  are  apparent  in  the  log-rms  coda  envelopes  in 
Figures  11  and  13.  The  coda  at  this  time  may  consist  of  weak, 
incoherent  shear  waves  produced  by  the  source  itself  and/or  from 


mode  conversion  or  forward  scattering  of  other  phases.  Like  Lg, 
shear  waves  may  be  more  sensitive  to  the  effects  of  laterally 
varying  structure  beneath  the  array  than  scattered  P  waves.  This 
is  because  S  waves  have  smaller  wavelengths,  and  hence,  may  be 
more  likely  to  be  scattered  by  small-scale  lateral  heteroge¬ 
neities  than  P  waves. 

2.5  Lg-TO-P  FORWARD  SCATTERING 

The  flat  part  of  the  codas,  observed  at  NORSAR  between  200 
and  340  seconds  after  P,  appears  to  be  produced  by  a  long  burst 
of  energy  that  is  not  associated  with  any  known  phase.  We  have 
indicated  in  Figure  13  the  expected  arrival  times  of  major 
phases,  and  the  flat-coda  does  not  appear  to  be  associated  with 
any  of  them.  The  major  phases,  PP  and  PcP,  arrive  well  ahead  of 
the  flat  part,  and  S  waves,  including  Sn,  would  be  expected  to 
arrive  at  NORSAR  after  the  time  indicated  for  S  in  Figure  13.  A 
third  possibility  is  the.  arrival  could  be  Pg,  although  to  our 
knowledge,  there  have  been  no  observations  of  Pg  at  this  great  a 
distance  (A=  38°).  From  close  analysis  of  these  codas  we  con¬ 
clude  that  the  flat  part  of  the  coda  begins  at  about  219  seconds 
after  P  and  ends  at  352  seconds  after  P.  These  times  correspond 
to  group  arrival  times  of  6.42  and  5.34  km/sec,  respectively, 
which  is  within  the  reported  range  of  Pg  velocities  of  5.2  to  6.0 
km/sec.  However,  as  pointed  out  by  (Pomeroy  et  al,  1982),  Pg  is 
usually  weak  or  unidentifiable  at  distances  beyond  5  to  8°  in 
stable  cratonic  continental  interiors,  such  as  eastern  United 
States  and  analogous  regions  in  Russia. 


In  an  earlier  study  (Baumgardt,  1984)  ,  we  suggested  that  Lg 
forward  scattering  to  P-waves  may  generate  many  of  the  pre-Lg  P- 
coda  waves.  If  we  hypothesize  that  the  flat  part  of  the  coda  is 


a  burst,  of  P-wave  energy  caused  by  Lg  scattering,  we  can  deter¬ 
mine  where  the  scattering  occurs  by  timing  the  burst  relative  to 
Lg  and  use  a  regional  travel-time  curve  to  find  the  distance  to 
the  scattering  point  based  on  the  P-to-Lg  time.  The  flat  part  of 
the  coda  appears  to  begin  and  end  roughly  at  549  and  416  seconds 
before  Lg,  respectively.  From  the  regional  travel  time  curve  of 
(Gupta  et  al,  1980),  shown  in  Figure  16,  we  see  that  P  to  Lg 
times  of  549  and  416  seconds  correspond  to  distances  of  23°  and 
29°  from  NORSAR  respectively.  In  Figure  17,  where  the  Lg  propa¬ 
gation  paths  from  Semipalatinsk  to  NORSAR  and  Graefenburg  are 
shown,  we  see  that  the  distance  range  of  23°  to  29°  from  NORSAR 
brackets  the  north-south  trending  Ural  Mountains. 

The  Lg-to-P  forward  scattering  mechanism  is  illustrated 
schematically  in  Figure  18.  The  Lg  wave  is  generally  viewed  as 
consisting  of  a  superposition  of  S  waves  multiply  reflected 
within  the  crust  with  near-critical  incidence  angles  at  the  Moho 
(Campilla  et  al,  1984) .  Thus,  for  plane-parallel  layering  in  the 
crust,  there  would  be  verly  little  leakage  of  energy  into  the 
mantle.  However,  if  the  crust-mantle  or  the  free-surface  inter¬ 
faces  are  dipping,  such  as  in  the  Ural  Mountain  region,  the 
incidence  angles  at  the  Moho  may  become  subcritical,  and  the  S 
waves  may  leak  energy  into  the  lower  mantle  by  mode  conversion  to 
P.  The  P-waves  would  then  arrive  as  precursors  to  Lg.  Also,  S 
waves  may  propagate  into  the  mantle  and  produce  an  S-wave  coda  at 
times  before  Lg  but  after  the  expected  arrival  of  S  from  the 
source.  Another  possibility  is  that  a  mountain  range  like  the 
Urals  may  have  a  number  of  deeply  penetrating  near-normal  faults 
which,  if  there  has  been  substantial  vertical  displacement  along 
the  faults,  may  produce  horizontal  impedance  contrasts  across  the 
fault  surfaces.  Thus,  the  horizontally  propagating  S  waves, 
which  compose  the  Lg-wave  packets,  may  impinge  upon  the  fault 
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surface  and  undergo  mode  conversion  to  P  as  a  result  of  this 
impedance  contrast. 

As  seen  in  Figure  17,  the  path  from  Semipalatinsk  to  the 
Graefenburg  just  misses  the  southern  tip  of  the  Urals  which 
explains  why  the  coda  flattening  is  not  observed  in  the 
Graefenburg  coda  plots  (Figures  8  and  9).  In  contrast  to  the 
NORSAR  codas,  the  Graefenburg  codas  decay  monitonically  with  time 
until  the  possible  arrival  of  Sn  waves  at  about  600  seconds  after 
P.  At  this  point,  the  coda  levels  begin  to  increase  and  then 
peak  after  the  arrival  of  Lg. 

We  conclude,  therefore,  that  the  flat  envelopes  from  about 
200  to  340  seconds  after  P,  recorded  at  NORSAR  for  Semipalatinsk 
explosions,  are  consistent,  on  the  basis  of  timing,  to  Lg  -  P 
forward  scattering  from  lateral  heterogeneities  in  the  Ural  Moun¬ 
tains.  The  evidence  is  still  circumstantial,  and  we  will  be 
analyzing  additional  data  for  events  on  either  side  of  the  Urals. 
However,  if  the  Lg-to-P  forward  scattering  hypothesis  holds  up, 
it  would  have  significant  implications  for  Lg  scattering  attenua¬ 
tion  determinations  at  teleseismic  distances  and,  perhaps,  at 
regional  distances  as  well.  Moreover,  the  low  standard  deviation 
of  the  flat  part  of  the  coda  indicates  that  measurements  of  coda 
magnitude  in  this  part  of  the  coda  may  be  more  stable  for  yield 
estimation  than  P-wave,  Lg  or  early  P-coda  measurements,  which 
all  have  higher  standard  deviations. 


3.0  RELATIVE  Lg  AND  P-CODA  MAGNITUDE  ANALYSIS 


3.1  INTRODUCTION 

In  this  section,  we  examine  the  conjecture  of  Sykes  and 
Cifuentes  (1983)  that  the  largest  Shagan  River  explosions  have 
nearly  the  same  yield.  Specifically,  we  focus  on  those  explo¬ 
sions  in  Table  1  which  have  body-wave  magnitudes,  as  determined 
by  Sykes  and  Cifuentes  (1983),  in  excess  of  6.0.  They  argue  that 
events  with  mb  greater  than  6.0  have  about  the  same  yield  and 
that  their  mbs  differ  because  of  biasing  effects  caused  by  dif¬ 
fering  geology  in  the  source  region.  Also,  it  is  conceivable 
that  differing  amounts  of  non-isotropic,  tectonic  component  for 
the  different  events  may  produce  variation  in  mb,  although  the 
extent  to  which  short-period  P-waves  are  affected  by  tectonic 
component  is  not  well  understood.  As  we  discussed  in  the  pre¬ 
vious  section,  P-coda  and  Lg-wave  measurements  appear  to  be  less 
sensitive  to  biasing  effects  than  P-wave  measurements.  Thus,  we 
expect  that  the  P-coda  and  Lg  amplitudes  for  the  largest  Shagan 
River  explosions  should  be  nearly  the  same  if  the  conjecture  of 
Sykes  and  Cifuentes  (1983)  holds  true. 

Log-rms  amplitudes  of  Lg  and  P-coda  waves  recorded  at  NORSAR 
and  Graefenburg  were  determined  using  equation  3  in  Section  2.0. 
P-coda  log-rms  amplitudes  were  averaged  over  ten  5  second  win¬ 
dows,  or  50  seconds,  starting  5  seconds  after  the  first-arrival  P 
time.  For  Lg,  log-rms  values  in  twelve  5  second  windows  or  120 
seconds,  were  averaged,  beginning  40  seconds  ahead  of  the  expec¬ 
ted  3.5  km/sec  group  arrival  time.  Since  we  are  only  interested 
in  relative  magnitude  estimates,  no  attenuation  correction  was 
made  for  either  the  Lg  or  P-coda  measurements  in  order  to  obtain 
absolute  magnitudes. 


3-1 


The  Lg  and  P-coda  measurements  for  NORSAR  are  discussed  in 
detail  in  Ringdal  (1983)  and  Baumgardt  (1983),  respectively.  In 
the  next  section  we  discuss  the  Graefenburg  measurements. 

3.2  GRAEFENBURG  Lg  AND  P-CODA  MEASUREMENTS 

Figure  3  shows  that  the  Graefenburg  array  has  three  3- 
component  sensors,  labeled  Al ,  B1  and  Cl.  Figure  19  shows  a  plot 
of  the  log-rms  amplitudes  of  Lg  on  the  transverse  component 
against  those  on  the  vertical  component.  Transverse  component 
traces  were  made  by  mathematically  rotating  the  north-south  and 
east-west  horizontal  components  into  the  back  azimuth  of  Semi- 
palatinsk  relative  to  Graefenburg.  Measurements  on  the  vertical 
and  transverse  components  from  all  three  sites  are  plotted  in 
Figure  19.  Each  trace  was  bandpass  filtered  between  0.6  and  3.0 
Hz  prior  to  making  measurements. 

Figure  19  shows  that  the  relative  log-rms  Lg  amplitude 
measurements  on  the  vertical  and  transverse  components  at  the 
three  sites  are  not  perfectly  consistent.  The  transverse  Lg  at 
the  Bl  site  is  between  0.2  to  0.4  log-rms  units  higher  relative 
to  the  vertical  Lg  as  compared  with  the  Al  and  Cl  sites.  The 
vertical  and  transverse  Lg  amplitudes  at  the  Al  and  Cl  sites  are 
more  consistent,  although  the  Al  amplitudes  are  somewhat  higher 
on  both  components  than  the  Cl  amplitudes.  Variations  in  geology 
between  the  sites  probably  accounts  for  the  differences  in  ampli¬ 
fication  of  transverse  relative  to  vertical  Lg .  Gupta  et  al 
(1982)  studied  variations  in  hor izontal-to-ver tical  Lg  amplitude 
ratios  in  the  eastern  United  States  and  found  that  the  ratio  is 
larger  at  soft-rock  sites  (H/Z  >2)  than  hard  rock  sites  (H/Z  <2). 
The  H/Z  ratios  at  the  Bl  three-component  site  are  as  large  as  1.6 
compared  with  average  H/Z  values  of  1  for  the  Al  and  Cl  sites. 
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We  thus  conclude  that  31  is  situated  on  softer  rock,  or  underlain 
by  more  layered  sediments,  than  A1  and  C.’  - 

We  have  also  noticed  that  the  noise  levels  on  the  B1  and  Cl 
instrument  are,  in  many  cases,  higher  than  those  recorded  by  the 
Al  instrument.  Because  of  this  fact  and  the  relatively  high 
correlation  in  the  vertical  and  transverse  amplitudes  at  Al,  we 
have  concentrated  on  measurements  on  the  Al  vertical  and  trans¬ 
verse  Lg  and  P-coda  for  the  relative  magnitude  analysis. 

Because  of  the  availability  of  broadband  data,  we  determined 

Lg  and  P-coda  log-rms  amplitudes  in  two  frequency  bands,  0.2  to 

1.0  Hz  and  0.6  -  3.0  Hz.  The  latter  covers  the  upper  tart  of  the 

mid-period  band  which  lies  between  the  traditional  low-  and  high- 

frequency  bands.  Lg  and  P-coda  measurements  were  made  in  both 

bands  on  both  the  vertical  and  transverse  components.  The 

i  i  i  i 

results  are  tabulated  in  Table  2,  where  mlgzr  migt  mcz'  and  mct 
refer  to  vertical-Lg,  transverse-Lg ,  ver tical-P-coda  and  trans- 
verse-P-coda  log-rms  amplitude  measurements,  respectively.  The 
superscript  i  =  l  and  2  corresponds  to  measurements  in  the  high- 
frequency  (0.6  to  3.0  Hz)  and  mid-frequency  (0.2  to  1.0  Hz) 
bands,  respectively. 

Figures  20  through  27  are  scatter  plots  of  each  of  the 
measurements  in  Table  2  against  the  m^  values  of  Sykes  and 
Cifuentes  (1983)  .  Also  shown  on  the  plots  are  the  results  of 
two-variable-in-error  regressions  with  A  =  1  and  A  =0  (Madansky, 
1959).  Each  observation  is  assumed  to  have  the  same  variance, 
and  A=  ffx/  a  y,  where  <jx  and  are  the  standard  deviations  of  the 
independent  and  dependent  variables,  respectively.  When  A=l,  we 
assume  equal  uncertainty  in  the  Graefenburg  and  Sykes  and  Cifu¬ 
entes  (1983)  measurements,  whereas,  A=0  assumes  that  all  the 
uncertainty  resides  in  the  Graefenburg  measurement  with  no  error 


FIGURE  20  SCATTER  PLOT  OF  VERTICAL,  HIGH-FREQUENCY 

Lg  LOG -RMS  AMPLITUDES  RECORDED  AT 

GRAEFENBURG  A1  VERSUS  mu 
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FIGURE  21  SCATTER  PLOT  OF  VERTICAL  MID- FREQUENCY 
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are  not  strongly  affected  by  tectonic  component.  Also,  the 
scatter  in  the  points  plotted  in  these  figures  does  not  appear  to 
be  caused  by  tectonic  component.  For  example,  in  Figure  20, 
events  11,  12,  and  14  fall  well  off  the  best  fitting  lines,  but, 
from  Table  3,  they  have  relatively  low  tectonic  components.  On 
the  other  hand,  events  4,  5,  7,  and  10,  which  have  large  tectonic 
components,  lie  closer  to  one  or  both  of  the  best  fitting  lines. 

Another  important  observation  which  can  be  drawn  from  these 
plots  is  that  the  slopes  of  the  best  fitting  lines  for  the  Lg 
measurements  are  quite  small,  ranging  from  0.6  to  0.8.  Interest¬ 
ingly,  however,  the  P-coda  measurements  have  higher  slopes  of  0.9 
to  1.0.  Baumgardt  (1984)  has  argued  that  early  vertical  P-coda 
waves  are  produced  by  P-scattering  beneath  the  source,  receiver, 
and  along  the  path  between  source  and  receiver.  Thus,  the  ampli¬ 
tudes  of  P  and  early  coda  waves  should  be  more  correlated  whereas 
P-  and  Lg-wave  amplitudes  may  be  more  independent. 

One  possible  cause  of  the  low  slopes  is  that  the  Lg  ampli¬ 
tudes  of  the  smaller  events  may  be  more  contaminated  by  noise 
than  those  of  the  larger  events.  Although  this  is  undoubtably  a 
problem  for  events  with  mbs  less  than  5.7  or  5.8,  our  sample  of 
data  only  include  two  events  in  this  magnitude  range.  While  it 
is  true  that  undersampling  of  lower  magnitude  events  means  that 
the  slope  is  not  well  contained  by  the  smaller  events,  noise 
contamination  cannot  be  the  principal  cause  of  the  small  slopes 
because  most  of  the  events  produce  Lg  amplitudes  at  Graefenburg 
which  are  well  above  the  noise  level. 

The  principal  cause  of  the  small  slopes  is  that  many  of  the 
events  with  different  m^  have  close  to  the  same  Lg  and  P-coda 
amplitude.  This  is  particularly  apparent  for  the  largest  events, 
with  m^  greater  than  6.0,  where  there  is  a  noticeable  flattening 


of  the  scatter  plots  in  Figures  20  through  27.  Notice  that  the 
flattening  is  evident  in  the  P-coda  scatter  plots,  even  though 
their  slopes  are  larger  than  those  for  Lg.  Clipping  does  not 
appear  to  be  the  cause  of  the  flattening  because  the  observed  Lg 
and  coda  amplitudes  are  well  within  the  total  dynamic  range  of 
the  Graefenburg  recording  system. 

To  further  illustrate  this  observation,  we  have  replotted 
the  scatter  diagrams  in  the  form  of  histograms  in  Figures  28 
through  35.  The  histograms  of  each  Lg-  and  P-coda  measurement 
are  plotted  on  the  right  and  the  corresponding  mb  histogram  is 
plotted  to  the  left  in  each  figure.  Examination  of  each  mb 
histogram  reveals  two  peaks  centered  at  about  5.95  and  6.15. 
Sykes  and  Cifuentes  (1983)  also  pointed  out  these  two  peaks  in 
the  magnitude  distributions,  and  suggested  that  all  the  events 
which  cluster  within  the  distribution  centered  at  the  higher  mb 
have  the  same  yield.  (Actually,  their  higher  mb  distribution  was 
centered  on  6.2  and  included  events  with  magnitudes  from  6.1  to 
6.28.  However,  because  our  data  sample  is  smaller,  we  have 
included  an  event  with  magnitude  of  6.033  in  the  larger  magnitude 
distribution  which  Sykes  and  Cifuentes  (1983)  placed  in  the 
distribution  centered  at  about  5.98.  However,  we  do  not  believe 
that  this  will  invalidate  the  following  arguments.)  The  arrows 
indicate  the  approximate  spread  of  this  distribution  which  is  0.2 
mb  units.  The  corresponding  range  in  Lg  and  P-coda  log-rms 
amplitudes  (i.e,  for  the  same  events)  is  denoted  by  the  arrows  on 
the  left  histograms  in  each  of  Figures  27  through  34.  In  most 
cases,  the  spread  in  the  Lg  or  P-coda  amplitudes  is  less  than 
half  that  of  the  corresponding  mb  values. 

Thus,  we  conclude  that  the  small  slopes  in  the  scatter  dia¬ 
grams  are  a  result  of  the  Lg  and  P-coda  amplitudes  clustering 
into  distributions  with  tighter  variances  than  those  of  the 


FIGURE  28  HISTOGRAMS  OF  VERTICAL,  HIGH*-  FREQUENCY  Lg  LOG-RMS  AMPLITUDES 
RECORDED  AT  GRAEFENBURG  AND  M,  FOR  SHAGAN  RIVER  EVENTS 


FIGURE  29  HISTOGRAMS  OF  VERTICAL  MID- FREQUENCY  Lg  LOG-RMS  AMPLITUDES 
RECORDED  AT  GRAEFENBURG  AND  m,  FOR  SHACAN  RIVER  EVENTS 


FIGURE  30  HISTOGRAMS  OF  TRANSVERSE,  HIGH-FREQUENCY  Lg  LOG-RMS  AMPLITUDES 
RECORDED  AT  GRAEFENBURG  AND  m,  FOR  SHAGAN  RIVER  EVENTS 


FIGURE  31  HISTOGRAMS  OF  TRANSVERSE,  MID- FREQUENCY  Lg  LOG-RMS  AMPLITUDES 
RECORDED  AT  GRAEFENBURG  AND  m,  FOR  SHAGAN  RIVER  EVENTS 


FIGURE  32  HISTOGRAMS  OF  VERTICAL  HIGH-FREQUENCY  P-CODA  LOG -RMS 

AMPLITUDES  RECORDED  AT  GRAEFENBURG  AND  m,  FOR  SHAGAN  RIVER 
EVENTS  b 


FIGURE  34  HISTOGRAMS  OF  TRANSVERSE.  HIGH-FREQUENCY  P-CODA  LOG-RMS 

AMPLITUDES  RECORDED  AT  GRAEFENBURG  AND  m,  FOR  SHAGAN  RIVER 
EVENTS  D 


corresponding  m^  values.  This  means  that  the  Lg  and  P-coda 
amplitudes  are  closer  together  than  the  corresponding  m^ 
values.  We  now  interpret  this  result  in  terms  of  the  relative 
yields  of  the  largest  Shagan  River  explosions. 

3.3  RELATIVE  MAGNITUDES  AND  YIELDS  OF  THE  LARGEST  SHAGAN 

RIVER  EXPLOSIONS 

In  Table  4,  we  summarize  the  relative  magnitude  values  for 
all  Shagan  River  explosions  with  m^s  greater  than  6.0.  In  this 
table  we  compare  the  relative  values  of  several  magnitude  esti¬ 
mates,  including  NORSAR  Lg  and  P-coda  measurements  of  Ringdal 
(1983)  and  Baumgardt  (1984),  respectively,  the  Graefenburg 
results  discussed  above,  the  Lg  magnitudes  of  Nuttli  (1984) ,  and 
the  body-wave  magnitudes  of  NEIS  and  Sykes  and  Cifuentes  (1983) . 
From  this  table  we  make  the  following  observations. 

i.  The  relative  NORSAR  Lg  and  P-coda  magnitudes 
for  these  events  vary  by  ±0.04  to  ±0.05  loga¬ 
rithmic  units  as  compared  with  ±0.08  to  ±0.10 
for  P-wave  magnitudes. 

ii.  The  Graefenburg  Lg  and  P-coda  measurements 
vary  by  ±0.03  to  ±0.05  logarithmic  units  as 
compared  to  ±0.07  units  for  the  corresponding 
m^  values. 

iii.  The  variation  in  Lg  magnitudes,  estimated  by 
Nuttli  (1984)  is  ±0.07  as  compared  with  ±0.04 
for  the  corresponding  m^  values. 

Observations  (i)  and  (ii)  indicate  that  the  range  of  Lg  and 
P-coda  magnitudes  for  the  largest  explosions  is  half  that  of  the 
m^  range.  Statistical  F  tests  reveal  that  these  observations  are 
significant  at  and  above  the  95%  confidence  level.  Assuming  a 
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TABLE  4 

MEANS  AND  STANDARD  DEVIATIONS  FOR 
MAGNITUDE  MEASUREMENTS  OF 
LARGEST  SOVIET  EXPLOSIONS, 
1976-1982 


Method 

Mean 

Standard 

Deviation 

No.  of 
Events 

1.  Single-Channel 

Lg  NORSARI 

3.11 

0.04 

5 

Multi-Channel 

Lg  NORSARTT 

3.01 

0.05 

5 

NETS  mK  Corres- 

ponding  Events 

6.12 

0.13 

5 

Corrected 

Corresponding 

Events 

6.143 

0.092 

5 

2.  Multi-Station 

Lg  WWSSN 

6.01 

0.07 

4 

NEIS  m.  Corres- 

ponding  Events 

6.18 

0.10 

4 

Corrected  mb 
Corresponding 

Events 

6.167 

0.036 

4 

3.  Single-Channel 

P-coda  NORSAR^ 

1.97 

0.05 

6 

Multi-Channel 

P-coda  NORSAR++ 

1.96 

0.06 

6 

NEIS  mh  Corres- 

ponding  Events 

6.17 

0.10 

6 

Corrected 

Corresponding 

Events 

6.160 

0.082 

6 
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magnitude-yield  slope  of  1,  the  corresponding  yield  ranges  near 
±150  kt  are  +17  kt  for  Lg  and  P-cod*.  measurements  as  compared  to 
±35  kt  for  P-wave  magnitudes. 

The  Nuttli  (1984)  network  Lg  results  are  not  consistent  with 
the  single-array  results  in  that  he  gets  a  larger  magnitude  range 
(±0.07)  for  Lg  than  for  m^  (±0.04).  A  possible  explanation  for 
the  discrepancy  is  that  the  number  of  stations  used  by  Nuttli  to 
estimate  the  average  Lg  magnitude  was  not  the  same  from  event  to 
event.  The  number  of  stations  used  ranged  from  2  to  6  stations. 
This  large  variability  of  the  number  of  stations  may  have  resul¬ 
ted  in  biased  average  magnitudes,  particularly  in  the  cases  where 
the  number  of  stations  is  small.  A  better  approach  would  have 
been  to  apply  the  Least  Squares  Matrix  Factorization  (LSMF)  ,  or 
analysis  of  variance,  method  to  estimate  average  magnitudes,  as 
did  Roundout  Associates  (1984) .  Also,  it  should  be  noted  that 
Nuttli' s  Lg  magnitudes  are  determined  by  measuring  the  single¬ 
point  peak  to  trough  amplitude  on  the  maximum  Lg,  whereas  we  have 
measured  log-rms  amplitudes  in  long  windows.  Log-rms  measure¬ 
ments  may  be  more  stable  than  single-point  measurements. 

Table  5  summarizes  the  results  of  the  m^  and  Ms  estimates 
for  the  largest  Shagan  River  explosions  (m^ >  6.0)  given  by  Sykes 
and  Cifuentes  (1983) .  The  relative  Ms  values  for  all  events  with 
m^  greater  than  6.0  and  including  events  of  high  and  low  tectonic 
component  have  a  large  standard  deviation  (±0.15),  whether  or  not 
tectonic  release  corrections  are  applied.  The  m^  values  for  the 
corresponding  events  has  a  standard  deviation  of  better  than  half 
(±0.07)  that  for  the  Ms  values.  In  the  case  of  low  tectonic 
release  events,  with  or  without  tectonic  release  corrections,  the 


TABLE  5 

SURFACE  WAVE  RESULTS  FOR  THE 
LARGEST  S HAGAN  RIVER  EXPLOSIONS 
(SYKES  AND  CIFUENTES,  1983) 


Method 

Mean 

Standard 

Deviation 

No.  of 
Events 

1.  Multi-Station  M 

WWSSN  All  Events  - 
Uncorrected 

4.004 

0.15 

9 

Multi-Station  M 

WWSSN  All  Events  - 
Corrected4- 

4.190 

0.15 

7 

Corrected 

All  Events 

6.158 

0.072 

10 

2.  Multi-Station  M 

WWSSN  Low  TR 
(LQ/LR  0.4) 
Uncorrected 

4.097 

0.069 

5 

Multi-Station  Ms 
WWSSN  Luw  TR 
(LQ/LR  0.4) 

Correct  ?d+ 

4.246 

0.072 

5 

Corrected 

Corresponding 

Events 

6.118 

0.080 

5 

+ 

* 


Ms  corrected  by  Sykes  and  Cifuentes  (1983)  for  pure  thrust 
faulting . 

Station  corrections  derived  by  analysis  of  variance. 


Ms  and  mb  ranges  are  more  comparable  (±0.07  to  ±0.08)  although 
still  about  double  that  which  we  obtain  for  Lg  and  P-coda 


measurements . 


Sykes  and  Cifuentes  (1983)  based  their  conjecture,  that  the 
largest  Shagan  River  events  are  of  nearly  the  same  yield,  on  the 
relative  Mg  values  of  only  a  few  events  of  low  tectonic  release. 
They  concluded  that  these  explosions  are  within  10  to  20  kt  of 
each  other.  Our  results  are  consistent  with  a  somewhat  larger 
yield  range  of  30  to  35  kt,  although  we  included  events  of  large 
tectonic  release  (LQ/LR  >0.4).  However,  based  on  the  Lg  and  P- 
coda  measurements,  we  concur  with  the  overall  conclusion  of  Sykes 
and  Cifuentes  (1983)  that  the  yield  range  of  the  largest  Shagan 
River  explosions  is  smaller,  by  at  least  a  factor  of  two,  than 
the  yield  range  indicated  by  the  mh  estimates. 
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coda,  where  the  spatial  standard  deviation  in  log-rms  amplitude 
in  5  second  windows  drops  to  a  minimum  of  less  than  0.1  units. 

5.  The  flat  coda  envelopes  from  about  200  to  340  sec¬ 
onds  after  P,  recorded  at  NORSAR  for  Semipalatinsk  explosions, 
are  consistent,  on  the  basis  of  timing,  to  Lg  -  P  forward  scat¬ 
tering  from  lateral  heterogeneities  in  the  Ural  Mountains.  The 
evidence  is  still  circumstantial,  and  we  will  be  analyzing  addi¬ 
tional  data  for  events  on  either  side  of  the  Urals.  However,  if 
the  Lg-to-P  forward  scattering  hypothesis  holds  up,  it  would  have 
significant  implications  for  Lg-scatter ing  attenuation  determina¬ 
tions  at  teleseismic  distances  and,  perhaps,  at  regional  dis¬ 
tances  as  well. 

6.  Measurements  on  Graefenburg  data  of  vertical  and 
transverse  log-rms  amplitudes  in  the  P-coda  and  Lg  phase  in  two 
frequency  bands  were  found  to  correlate  well  with  network  mb. 
However,  the  slopes  of  the  regression  of  Lg  amplitude  versus  mb 
were  small  ranging  from  0.6  to  0.8.  The  slopes  of  the  P-coda 
regressions  were  higher,  ranging  from  0.9  to  over  1.0.  This 
supports  the  idea  that  early  P-coda  waves  are  generated  by  P 
scattering . 

7.  Lg  and  P-coda  log-rms  amplitudes  do  not  appear  to 
be  strongly  affected  by  the  tectonic  component  as  revealed  by 
surface-wave  studies. 

8.  Relative  P-coda  and  Lg-magnitude  analysis  of  the 
largest  Shagan  River  explosions,  with  mb  >  6.0,  using  NORSAR  and 
Graefenburg  data  supports  the  conjecture  of  Sykes  and  Cifuentes 
(1983)  that  these  explosions  have  nearly  the  same  yield. 


9.  Of  course,  we  cannot  address  the  other  contention 
of  Sykes  and  Cifuentes  (1983)  that  these  explosions  do  not  exceed 
the  150  kt  threshold,  since  we  did  not  estimate  absolute  Lg  and 
P-coda  magnitudes  which  are  calibrated  against  the  yields  of 
events  with  known  yield.  Nuttli  (1984)  has  done  this  and  has 
obtained  higher  yields,  some  in  excess  of  200  kt,  for  the  largest 
explosion.  However,  we  question  the  precision  of  his  estimates 
because  the  standard  deviation  between  events  of  his  network- 
averaged  Lg  magnitudes  for  the  largest  Shagan  River  explosion  do 
not  agree  with  our  single-site  estimates  at  both  NORSAR  and 
Graefenburg.  Perhaps  his  magnitudes  have  bias  because  of  the 
variable  number  of  stations  used  to  estimate  average  magnitudes 
and  his  use  of  single-point  rather  than  rms  amplitudes  to 
determine  Lg  magnitudes. 

10.  Single-station  Lg  and  P-coda  measurements  appear  to 
be  significantly  more  precise  for  relative  yield  estimation  than 
any  other  method.  If,  in  fact,  the  largest  Shagan  River  explo¬ 
sions  have  about  the  same  yield,  then  we  can  conclude  that 
single-station  P-coda  and  Lg  magnitudes  have  a  precision  of  about 
±35  kt,  or  a  factor  of  1.6,  for  relative  yield  estimation  com¬ 
pared  with  network  m^  precision  of  ±60  kt  or  a  factor  of  2.3. 

These  conclusions  are  based  on  a  very  small  data  sample  and, 
therefore,  must  still  be  considered  preliminary.  We  are  current¬ 
ly  collecting  more  data  from  both  arrays  for  more  recent  events 
which  should  give  us  a  more  statistically  meaningful  sample  for 
relative-magnitude  analysis.  Also,  we  will  be  collecting  data 
from  other  digital  stations,  such  as  the  SRO  and  DWWSSN  stations 
in  Europe  and  Scandinavia.  We  plan  to  combine  this  data  with  the 
array  data  to  investigate  the  precision  of  average  Lg  and  P-coda 


We  recommend  that  more  attention  be  given  to  broadband 
recordings  of  Lg,  such  as  from  Graefenburg,  intermediate  band 
data  from  DWWSSN  and  RSTN  stations  and  from  the  newly  installed 
NORESS  broadband  3-component  sensor.  We  find  that  Lg  excitation 
is  stronger  in  the  intermediate  band.  However,  noise  is  often 
correspondingly  higher  in  this  band  as  well,  so  methods  for  noise 
correction  need  to  be  investigated. 

Finally,  Lg  measurements  at  different  frequencies  need  to  be 
compared  to  Ms  measurements,  with  and  without  tectonic  release 
corrections.  Our  preliminary  examination  of  tectonic  release 
corrections  given  by  Sykes  and  Cifuentes  (1983)  suggest  that  they 
are  not  very  effective  for  improving  yield  estimation  precision, 
assuming  the  largest  Shagan  River  explosions  are  of  nearly  the 
same  yield.  It  should  be  remembered  that  the  Sykes  and  Cifuentes 
(1983)  conclusions  were  based  almost  entirely  on  Ms  measurements 
for  events  with  very  low  tectonic  component,  and  they  warn  a- 
gainst  using  Ms  for  yield  estimation  of  events  whose  LQ/LR  ratios 
much  exceed  0.4.  More  recent  Ms  corrections,  based  on  moment 
tensor  inversions,  should  be  compared  with  Lg  measurements. 
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